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ABSTRACT 



We present measurements of ratios of elements of the scattering matrix of 
Martian analogue palagonite particles for scattering angles ranging from 3° to 
174° and a wavelength of 632.8 nm. To facilitate the use of these measurements 
in radiative transfer calculations we have devised a method that enables us to 
obtain, from these measurements, a normalized synthetic scattering matrix 
covering the complete scattering angle range from 0° to 180°. Our method 
is based on employing the coefficients of the expansions of scattering matrix 
elements into generalized spherical functions. The synthetic scattering matrix 
elements and/or the expansion coefficients obtained in this way, can be used to 
include multiple scattering by these irregularly shaped particles in (polarized) 
radiative transfer calculations, such as calculations of sunlight that is scattered 
in the dusty Martian atmosphere. 



Keywords: scattering ; dust ; Martian analogue palagonite ; expansion coef- 
ficients ; polarization 
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1 Introduction 



Dust from the Martian surface is regularly swept up by winds to form lo- 
cal, regional, or sometimes even planet-wide dust storms. The airborne dust 
particles scatter and absorb solar radiation and are therefore very important 
for the thermal structure of the thin M artian atmosphere and for the tem- 



perature of the Martian surface (see e.g. iBougher et all 120061 ; iGurwell et al. 



20051 ; ISmithl . 12004 ; iGierasch and Goodyl . 1 1 9 721 . and references therein). The 



interaction between radiation and the dust particles thus has to be taken into 
account when studying, for example, the local and global climate on Mars. In 
particular, one has to account for the spatial distribution of the dust particles, 
their number density, and their optical properties. For a given wavelength, the 
optical properties of the dust particles depend on their composition, sizes and 
shapes. 



Despite the important role of dust particles in the Martian atmosphere, sur- 
prisingly little i s know n about their optical properties (for an overview, see 



Korablev et all 120051 ) . Consequently, in radiative transfer calculations that 
are used to interpret space or ground-based observations of Mars, various 
assumptions are made regarding the dust optical properties. In particular, 
it is common to assume homogeneous, s pheric al or spheroidal dust parti- 



cles ( Pollack et all Il995l ; iTomasko et al. 



19991 ) . although dust particles on 



Earth are known to be irregularly shaped. The optical properties of homo- 
geneous, sp herical particles can straightforwardly be calculated us ing Lorenz- 
Mie theory (Ivan de Hulstl . 119571 ; Ide Rooij and van der Stapl . Il984 ) , and those 
of homogeneous, spheroidal pa r ticles using e.g. the so -called T-matrix method 



( IMishchenko and Travisl . 1 19941 ; iDubovik et all 120061 ). The optical properties 



of spherical particles can, however, differ significantly from those of irregularly 
shaped particles, even if their composition and/or size distribution is the same. 
Therefore, assuming spherical instead of irregularly shaped particles in radia- 
tive transfer calculations that are used for example to analyze observations, 
can lead to significant errors in retrieved atmospheric parameters, such as the 
dust optical thickness and /or dust particle size distributions (for a discu ssion 
on such errors, see e.g. iDlugach and Petroval . 120031 ; iDlugach et all 120021 ). 



For irregularly shaped particles, which are very common in nature, the scat- 
tering matrix elements can in principle be calculated with numerical methods 
such as those b ased on the so-called Dis crete Dipole Approximation (DDA) 
method (see e.g. lDraine and Flat aul . 1 19941 ) . However, the vast amounts of com- 
puting time that such numerical calculations require make it at least very 
impractical to calculate complete scattering matrices for a sample of par- 
ticles with various (irregular) shapes and various sizes, in particular if the 
particles are large compared to the wavelength of the scattered light. Alterna- 
tively, employing geometrical optics for unrealistically spiky particles com- 
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bined with an imaginary part of the refractive index that is rather small 
compared to the typical values used in the literature, appears to be use- 
ful to re produce the scattering b ehavior of irregularl y shaped mineral par - 



ticles (see Nousiainen et all, 2003h. As suggested by e. g. Tomasko et al.l (jl999l ) 



Dlugach and Petrova fl2003l ) : IWolfT and Clancvi (l2003h . a more practical method 



to obtain elements of the scattering matrix for an ensemble of irregularly 
shaped particles is to measure the elements in a laboratory. Note that mea- 
sured scattering matrix elements are also essential for validating numerical 
methods and approximations. 



In this article, we present measurements of ratios of elements of the scattering 
matrix of irregularly shap ed, randomly orien ted Martian analogue palago- 
nite particles, described by iBanin et al.l (119971 ). as functions of the scattering 
angle. The material palagonite is believed to be a reasonable, but not per- 
fect, analogue for the Martian surface and atmospheric dust particles. Terres- 
trial palagonite particles (i.e. terrestrial weathering produc ts of basaltic ash or 
glass ) have been put forwar d as Martian dust analogues (jEvans and Adamsl . 
198ll : lRoush and Belj[l995h because of spectral similarities observed with vis- 



ible and near-infrared spectroscopic observati ons using both Earth-based tele - 
scopes and several Mars orbiting spacecraft (ISingerl . Il985l ; ISoderbloml . Il992l ). 



The measurements have been performed using a HeNe laser, which has a 
wavelength of 632.8 nm, and for scattering angles, G, ranging from 3° (near- 
forward scattering) to 174° (near-backward scattering). Other examples of 
such measurements for irregularly shaped mineral pa rticles obtain e d wit h the 
same experimental set-up have been reported by e.g. Volten et al.l ( 2005 ) and 
references therein. Since we have measured ratios of all (non-zero) elements 
of the scattering matrix as functions of the scattering angle, our results can 
be used for radiative transfer calcul ations that include multiple scattering and 
polarization. As described by e.g. Lacis et aD ( 1998 ). ignoring polarization, 
i.e. using only scattering matrix element (1,1) (the so-called phase function), 
in multiple scattering calculations, induces errors in calculated fluxes. The use 
of only the phase function should be limited to single scattering calculations 
for unpolarized incident light. 



A practical limitation of our experimental method is that we cannot measure 
close (< 6°) to the exact backscattering direction (0 = 180°), because there 
our detector would interfere with the incoming beam of light, nor close (< 3°) 
to the exact forward scattering direction (0 = 0°), because there our de- 
tector would intercept the unscattered part of the incident beam. These two 
scattering directions are, however, important for radiative transfer applica- 
tions. This holds in particular for the near-forward scattering direction, since 
a significant fraction of the light that is incident on a particle is generally scat- 
tered in the near-forward direction. A solution to the lack of measurements 
in the near-forward and near-backward scattering directions is to add arti- 
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ficial data points. For the near-forward scattering direction, where a strong 
peak in the phase function is expected, one can add artificial data points 
calculat ed using e.g. Lorenz-Mi e calcu l ations. Th i s has been done before, 
e.g. b y iKahnert and Nqusiainenl (120071 ); iLiu et al.l (120031 ); IVeihelmann et al. 



(l2004h: iHerman et all fl2005h . In this article, we use a method similar to 



that of ILiu et al.l (120031 ). that was also used by e.g. iMunoz et al.l (120071 ). but 



we extend this by using expansion coefficients which result from a Singu- 



spherical functions ( 


Gel'fand et al.. 


1963 




Hovenier and van der Mee 




1983: 


de Rooij and van der Stapl. 


1984; 


Hovenier et al.. 


2004). With the added ar- 



tificial data points and the expansion coefficients, we construct a so-called 
synthetic scattering matrix, which is normalized so that the average of the 
synthetic phase function over all directions equals unity and covers the whole 
scattering angle range (i.e. from 0° to 180°). 

Tables of the measurements, the synthetic scattering matrix elements and the 
expansion coeffi cients will be available fro m the Amsterdam Light Scattering 
Databas^f]. See Volten et al. ( 2005 . 2006 ) for a description of this database. 



The structure of this article is as follows. In Sect. El we describe the microphys- 
ical properties of our Martian analogue palagonite dust particles. In Sect. El 
we define the scattering matrix, describe the experimental set-up, and present 
the measurements and an auxiliary scattering matrix. In Sect. HI we intro- 
duce the expansion coefficients, describe the Singular Value Decomposition 
fitting method, and present the derived expansion coefficients and synthetic 
scattering matrix. In Sect. [51 finally, we summarize and discuss our results. 



2 Martian analogue palagonite particles 



Palagonite is a fine-grained weathering product of basaltic glass. At visible 
wavelengths it has a refractive index, m, typical for silicate materials, i.e . 
Re(m ) is about 1.5 and Im(m) is in the range 10 -3 to 10~ 4 ( Clancy et all 



19951 ) . Palagonite contai ns a considerab l e frac tion (about 10% by mass) of 



iron (III) oxide (Fe 2 3 ) (IRoush and Belli . Il995l ). which gives Mars its reddish 
color. 



The sa mple in this study is sample 91-16 that is described in detail by lBanin et al 
(119971 ). Note that there is another Martian analogu e palagonite sample de - 



scribed in the literature, namely sample 91-1 (see IRoush and Belli . Il995l ) 



Palagonite sample 91-1 appears to contain more sodium than sample 91-16 



2 The Amsterdam Light Scattering Database 
http:/ /www. astro. uva.nl/scatter 



is 



located at: 
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because of evaporation and deposition of salt due to the proximity of its re- 
trieval site to the Pacific Ocean. Palagonite sample 91-16 was retrieved at the 
top of Hawaii's Mauna Kea volcano, about 4 km above sea level, where it was 
formed in a semi-arid environment likely associated with ephemeral melting 
water from ice. Hence, sample 91-16 is considered to be the better alternative 
for Martian dust of the two Martian palagonite analogues. 



Before using sample 91-16 in our light scattering experiment, we removed the 
millimeter-sized particles by using a sieve with a 200-yum grid width, to avoid 
clogging the aerosol generator. Figure [T] shows an image of the palagonite par- 
ticles obtained with a scanning electron microscope (SEM). This image clearly 
shows the irregular shapes of the palagonite particles. It should be noted that 
SEM images are not necessarily representative of the size distribution of the 
particles. The normalized projected-surface-area distribution of the dust par- 
ticles was measured by using a laser particle sizer that is based on diffraction 
without maki ng assumptions about the r efract ive indices of the materials of 
the particles (IKonert and Vandenberghd . 119971 ). From the projected-surface- 
area distribution, we derive the number distribution and the volume distribu- 
tion of the particles because these distributions are often required for numerical 
applications. Figure [2] shows the normalized number, volume, and projected- 
surface-area distributions of our Martian analogue palagonite particles as func- 
tions of logr, with r the radius of a projected-surface-ar ea equivalent spher e 
(for details on these size distributions, see Appendix A of lVolten et al.L 120051 ) . 
The number distribution of our palagonite particles was approximated by a 
log- normal distribution, yielding an effective radius, r eS , of 4.46 jttm and an 
effective variance, v e s, of 7.29. Note th at is a dimensionless p arameter. For 
precise definitions of r e s and v e s see lHansen and Travisl (119741 ). Eqs. (2.53) 
and (2.54), respectively. 



We are well aware of the fact that the sizes of real Martian dust particles can 
be very different from those in our sample. Indeed, sizes of dust particles on 
Mars will probably vary from location to location, and from time to time, 
especially when in local or global storms, dust particles are lifted up from 
the surface to be deposited somewhere else: depending on the atmospheric 
turbulence, the particles in the Martian atmosphere could have very different 
size distributions than those on the surface. 



The effective radius of 4.46 /xm of our sample particles is a factor of 2 to 3 
larger than the va l ues p ut forward for the effective radius of Martian dust by 
Wolff and Clancyl ( 2003 ). who analyzed observations by the The rmal Emission 



Spectr ometer (TES) on-board the Mars Global Surveyor, and by lPollack et al 



(119951). who analyz e d obs ervations performed by the Viking Lander. In partic- 
ular, iPollack et al.l (119951 ) derived an effective r adius of 1.85 ± 0.3 /iin. From 
Pathfinder me asurements. iTomasko et al.l (119991 ) derived an effective radius of 
1.6 ± 0.15 /iin. iLemmon et al.l (120041 ) derived values from observations by the 
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Mars Explora t ion R overs Spirit and Opportun ity that are similar to those of 
Pollack et all (Il995h and lTomasko et all fjl999h . 



Although our sample particles thus seem to be rather large, it should be 
noted that particle sizes as derived from observations will depend on the 
observing method, e.g. looking at diffuse skylight or at the surface, as well 
as on the retrieva l meth od. In particular, according to numerical similations 
by iDlugach et al.l (120021 ). effective radii that are derived for spheroidal dust 
particles at visibile wavelengths, under the assumption that these particles 
are spheri c al, ca n be significantly undere stimated. At infrared wavelengths, 



Min et al.l (120031 ) and iFabian et al.l (120011 ) show that absorption and emission 
processes, even in the small size parameter regime (i.e. 27rr e g/A < 1), depend 
on the particle shape, too. Clearly, because Martian dust is expected to show 
a great variety in microphysical properties, our results should simply be re- 
garded as an example of what can be expected for the scattering properties 
of irregularly shaped particles. 



3 The scattering matrix 



3. 1 Definition of the scattering matrix 



The flux and state of polarization of a quasi-monochromatic beam of light can 
be described by means of a so-called flux vector. If such a beam of light is scat- 
tered by an ensemble of randomly oriented particles, separated by distances 
larger than their linear dimensions and in the absence of multiple scattering 
as in our experimental set-up (see Sect. 13. 2p . the flux vectors of the incident 
beam, 7r$ (A), and scattered beam, 7r<fr(A, 9), are for each sc a ttering direc 



tion, related by a 4 x 4 matrix, as follows (Ivan de Hulst 
20061 ): 



1957; Volten et al. 



*(A,e) 



A 2 

47T 2 L> 2 





F\2 


F13 


F»\ 


Fn 


F22 


F23 


F24 


—F13 


—F23 


F33 


F34 


F u 


F24, 


—F34 


Fuj 



*o(A), 



where the first elements of the column vectors are fluxes divided by ir and the 
other elements describe the state of polarization of the beams by means of 
Stokes parameters. Furthermore, A is the wavelength, and D is the distance 
between the ensemble of particles and the detector. The scattering plane, i.e. 
the plane containing the directions of the incident and scattered beams, is the 
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plane of reference for the flux vectors. The matrix, F, with elements Fij is 
called the scattering matrix of the ensemble. The scattering matrix elements 
are dimensionless, and depend on the number of the particles and on their 
microphysical properties (size, shape and refractive index), the wavelength 
of the light, and the scattering direction. For randomly oriented particles, 
the scattering direction is fully described by the scattering angle 0, the angle 
between the directions of propagation of the incident and the scattered beams. 

According to Eq. (CQ), a scattering matrix has in general 10 different matrix 
elements. For randomly oriented particles with equal amounts of particles 
and their mirror particles, as we can assume applies for the particles of our 
ensemble, the four elements F 13 (0), Fu(Q), -F 23 (8 ) , an d -^24(6)) are zero over 
the entire scattering angle range (see Ivan de Hulstl . 119571 ). This leaves us only 
six non-zero scattering matrix elements, as follows 



F(0) 



F n {G) F 12 (0) 
F 12 (Q) F 22 (Q) 
F 33 (9) F 34 (9) 
-F 34 (6) F 44 (9) 



(2) 



where |i^(e)/F u (e)| < 1 flHovenier et all Il986h . 



For unpolarized incident light, matrix element Fn(O) is proportional to the 
flux of the singly scattered light and is also called the phase function. Also, 
for unpolarized incident light, the ratio —Fi 2 (Q) / Fu(Q) equals the degree 
of linear polarization of the scattered light. The sign indicates the direction 
of polarization: a negative degree of polarization indicates that the scattered 
light is polarized parallel to the reference plane, whereas a positive degree of 
polarization indicates that the light is polarized perpendicular to the reference 
plane. In calculations for fluxes only and where light is scattered only once, 
i*ii(0) is the only matrix element that is required. Ignoring the other matrix 
elements, and hence the state of polarization of the light, in multip le scattering 



calcu l ations, usually leads to errors in calcula t ed flu xes (see e.g. lLacis et al. 



19981 : Moreno et all 12002 : IStam and Hovenierl . 120051 ) 



3.2 The experimental set-up 



Our measurements have been performed with the ligh t scattering e xperi 



ment located in Amsterdam, the Netherlands (see e.g . IVolten et all 2005 



Hovenier et all 120031 : iHovenierl . l2000l : IVolten et all l200ll : iMunoz et all l2000h . 



Figure [3] shows a sketch of the experimental set-up. In our experimental ap- 
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paratus, we use a HeNe laser (A=632.8 nm, 5 mW) as a light source. The laser 
light passes through a polarizer and an electro-optic modulator. The mod- 
ulated light is subsequently scattered by an ensemble of randomly oriented 
particles from the sample, located in a jet stream produced by an aerosol 
generator. The scattered light may pass through a quarter-wave plate and 
an analyzer, depending on the s cattering matrix element of interest (for de- 



tails see e.g. IVolten et al.l . 120051 ). and is then detected by a photomultiplier 
tube which moves in steps along a ring with radius D (see Eq. (JTJ) around 
the ensemble of particles; in this way a range of scattering angles from 3° 
(nearly forward scattering) to 174° (nearly backward scattering) is covered in 
the measurements. 

We cannot measure close (< 3°) to the exact forward scattering direction, 
because there our detector would intercept the unscattered part of the incident 
beam, nor can we measure close (< 6°) to the exact backscattering direction, 
because there our detector would interfere with the incoming beam of light. 

A photomultiplier placed at a fixed position (i.e. at a fixed scattering angle) 
is used to correct the measured scattered fluxes for time fluctuations in the 
particle stream. It can safely be assumed th at during the measur ements, the 



particles are in the single scattering regime (jHovenier et al.l . 120031 ) 



Due to the lack of measurements between 0° and 3° and between 174° and 180°, 
we cannot measure the absolute angular dependency of the phase function, 
e.g. normalized to unity when averaged over all scattering directions. Instead, 
we normalize the measured phase function to unity at a scattering angle of 
30°. We present the other scattering matrix elements divided by the original 
measured phase function. We thus present ratios of elements of the scattering 
matrix instead of the elements themselves. 



3.3 Measurements 



Figure H] shows the six measured not identically zero (cf. Eq. [2]) ratios of 
elements of the scattering matrix of the Martian analogue palagonite parti- 
cles as functions of the scattering angle B, together with the experimental 
errors. We have verified that the measured ratios of the elements of the scat- 



tering matrix satisfy the Cloude coherency matrix test (jHovenier et al.l . 120041 ) 
within the experimental errors. And we verified that the other measured ratios 
of the elements of the scattering matrix, i.e. F 13 (0)/Fn(0), F 2 3(Q) / F U (Q) , 
-Fi 4 (9)/Fn(6), and i 7 24(©)/-^n(@), do not differ from zero by more than the 
experimental errors (see Eq. ([2])). 

To illustrate the influence of the particle shape on the scattering behavior 
of the palagonite particles, the measurements in Fig. @] are presented together 
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with results of Lorenz-Mie calculations (Ivan de HulstL 1 1 95 7t I de Rooij and van der Stap 
19841 ) for homogeneous, optically nonactive, spherical particles at a wavelength 
of 632.8 nm. For the Lorenz-Mie calculations we employed the number size 
distribution, n(r), derived from the measured projected-surface-area distribu- 
tion, and the refractive index was fixed to m — 1.5 + 0.0005i (cf. Sect. [2] and 

Fig. E]). 

As can be seen in Fig. HI the measured phase function, i.e. F 11 (0)/F 11 (3O°), 
of the irregularly shaped Martian analogue palagonite particles covers almost 
three orders of magnitude between = 3° and O = 174°, with a strong peak 
towards the smallest scattering angles (the so-called forward scattering peak) 
and a smooth drop-off towards the largest scattering angles. The measured 
phase function is very flat for scattering over intermediate (70° < G < 150°) 
and large (6 > 150°) scattering angles. The relatively flat appearance of the 
phase function of the palagonite particles at large scattering angles appears 
to be a general behavior for (terrestrial) ir regularly shaped m i neral particles 



with moder ate refractive indices (see e.g. IVolten et all 1200 ll ; iMunoz et al. 



2000l . l200ll ). Our palagonite phas e function res e mbles the phase functions 



measu red in-situ with the Viking (jPollack et all Il995l ) and Pathfinder mis- 



sions ( ITomasko et allll999l ). In Sect. El we make a more detailed compar ison 



between our phase function and those presented by lTomasko et al.l (119991 ) . 



As mentioned before (see Sect. 13. ip . the ratio — Fi 2 (0)/-Fn(9) represents the 
degree of linear polarization of the singly scattered light for incident unpolar- 
ized light. For the irregularly shaped palagonite particles, Fig. H] shows that 
this ratio has a characteristic (positive) bell shape at intermediate scatter- 
ing angles and a small negative branch for > 160°. For scattering an- 
gles larger than about 140°, the scattering angle dependence of our mea- 
sured ratio — F 12 (0)/-F 11 (0) resembles Earth-based observations of the plan- 
etary phase angle dependence of the degree of lin ear polarization of Mars 
( Dlugach and Petrova . 20031 ; Shkuratov et all 2005 ) (the planetary phase an- 
gle equals 180° — for single scattering). This suggests that the polariza- 
tion opposition effect that is observed at small phase angle s for most solid 



solar system bodies (see e.g. iRosenbush and Kiselevl . 120051 . and references 



therein) can be explained, at least partly, by single scattering by small ir- 
regular particles. Her e , it sh ould be noted that the obse rvations discussed by 
Dlugach and Petroval (120031 ) and IShkuratov et al.l (120051 ) pertain to light that 
has been scattered in the Martian atmosphere combined with light that has 
been reflected by the surface. It is thus not purely representative for airborne 
dust particles. 



The most striking difference between the measured and the calculated ratios 
— i 7 'i2(0)/i 7 ii(0) (Fig.HJ) is their sign, hence the direction of polarization of the 
scattered light for unpolarized incident light. The irregularly shaped particles 
mostly yield scattered light polarized perpendicular to the reference plane, 
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while the spherical particles yield scattered light polarized parallel to this 
plane. Another difference is that for the irregularly shaped particles, ratio 
— F 12 (0)/i 7 i 1 (0) is a smooth, almost featureless function of 6, while for the 
spherical particles, the ratio shows strong angular features, especially at large 
scattering angles. 



Scattering matrix element ratio i^22(©)/-Pii(@) is often used as a measure for 
the non-sphericity of the scattering particles, since for homogenous, optically 
inactive spheres, this ratio equals unity at all scattering angles. As can be seen 
in Fig. HI for the irregularly shaped palagonite particles, F 2 2(Q)/Fu(Q) devi- 
ates significantly from unity at all but the smallest scattering angles. Indeed, 
with increasing scattering angle, it decreases to slightly below 0.4 at G ~ 130°, 
and then increases again to 0.5 when approaches 180°. The scattering angle 
dependence measured for the palagonite particles is sim ilar in shape t o tha t 
reported for irregularly shaped mineral aeros ol particles (IVolten et all l200ll ). 



an d for e.g. var i ous ty pes of volcanic ashes (IMunoz et all 120041 ). According 



to IVolten et al.l (120011 ). the minimum value at intermediate scattering angles 
and the maximum value at the largest scattering angles are affected by the 
size and refractive index of the particles. 



Another indication for the shape of the scattering particles are the ratios 
F 33 (<d)/F u (<d) and F 44 (0)/F n (0). As can also be s een in Fig. H for ho - 
mogenous, optically inactive spheres, F 33 (Q) = F 44 (0) (IHovenier et all 120041 ) . 
whereas we find significant differences between the measured -F 44 (0)/F 1:L (0) 
and F 33 (Q) / F n (Q) for the palagonite sample. The ratio F 33 (Q)/Fn(Q) is zero 
at a smaller scattering angle than _F 44 (0)/.Fii(0), and has a lower minimum (- 
0.5 versus -0.2). Indeed, for the irregularly shaped particle s, these ratios show 
an ap parently typical behavior for non- spherical particles ( Mishchenko et all 
20001 ). namely, at large scattering angles, F 44 (0) / 'Fn(G) is larger than F 33 (Q) / F U (Q) . 



Finally, scattering matrix element ratio F 34 (0) / F n (0) of the irregularly shaped 
particles shows a shallow bell shape with slightly negative branches for < 30° 
and for > 165°. This scattering angl e dependence is commonly found for ir- 



regul arly shaped silicate particles (e.g. IVolten et all 120051 . 1200 ll ; IMunoz et al.l 



20001 ). Interestingly, whereas for the irregularly shaped particles, F 34 (0)/Fn(0) 
is very similar to — F 12 (0)/-Fn(0), for the spherical particles, these ratios dif- 
fer strongly from each other, both in sign and in shape, as can be seen in 
Fig. 1 



Comparison between the measured and the calculated scattering matrix ele- 
ment ratios in Fig. H] supports the idea that scattering by non-spherical par- 
ticles generally leads to smoother functions of the scattering angle than scat- 
tering by spherical particles. This smooth scattering behavior by irregularly 
shaped particles proves to be very difficult to simulate nu merically without 



taking into account the irregular shape of the particles (jNousiainen et al. 
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20031 ; iNousiainen and Vermeulenl . 120031 ; iKokhanovskyl . 120031 ) . 



An electronic table of the measured ratios of the elements of the scatter- 
i ng matrix wi l l be availab le from the Amsterdam Light Scattering Database 
flVolten et all 120051 . l2006h . 



3.4 The auxiliary scattering matrix 



It appears to be difficult to directly use the measured ratios of elements of 
the scattering matrix in radiative transfer calculations, because of the lack of 
measurements below = 3° and above G = 174°. In particular, it would be 
interesting to have the forward scattering peak in the phase function since 
it contains a large fraction of the scattered energy (see Fig. @J, and is thus 
very important for the accurate modelling of scattered light in e.g. planetary 
atmospheres. In addition, the lack of measurements at small and large scatter- 
ing angles inhibits the normalization of scattering matrix elements such that 
the average of the phase function over all scattering directions equals unity. 
With such a normalization and a value for the single scattering albedo of the 
scattering particles, one could model the absolute amount of radiation that is 
scattered in a given direction. 

To facilitate the use of the measured ratios of elements of the scattering matrix 
in radiative transfer calculations, we construct from these a so-called auxiliary 
scattering matrix, F au , for which holds (for i,j = 1 to 4 with the exception of 

i = J = l), 

^7( @ ) = J^§) F n (3) 



where the auxiliary phase function is equal to 

(©) = jr^fn ( 30 °) for 3°<0<174°. (4) 



This auxiliary phase function is normalized according to 

i- 1^(0)^=1, (5) 



■Lit 



where du is an element of solid angle. Combining Eq. (00) and Eq. (jSJ) and 
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setting ^(30°) equal to l/C leads to 



-I 

4tt J 

4-7T 



Fn(30°) 



du = C. 



(6) 



where C is a normalization constant. This constant can in principle be ob- 
tained by evaluating the integral on the left-hand side of Eq. ([6]), provided 
the function to be integrated is known over the full range of scattering angles. 
Therefore, we added artificial datapoints at and 180 degrees to the measured 
values of Fn(6)/Fn(30 o ). At = 180°, the smoothness of the measured phase 
function allows us to si mply add an artifi cial data point to Fn(6)/Fn(30 o ) 
by spline extrapolation (IPress et all Il992l ) of the measured data points. 



Adding an artificial data point to the measured F 11 (0)/F 11 (3O°) at 6 = 0°, 
is more complicated. Numerical tests with the calculated phase function for 
the hypothetical homogeneous, spherical palagonite particles (see Fig. U|) show 
that extrapolation of the calculated p hase function at < 3° towards 6 = 0°, 



using e.g. splines (IPress et all Il992l ). fails to reproduce the strength of the 



calculated forward scattering peak. We thus decided not to extrapolate the 
measured phase function from = 3° towards = 0°. Instead we add an ar- 
tificial data point to the measured Fn(0)/Fn(3O°) at = 0° using the phase 
function that we calculated for the projected-surface-area equivalent, homo- 
geneous, spheri cal particl e s. Th e rati onale for this appro ach, which is similar 
to that used by lLiu et al.l (120031 ) and lMunoz et al.l (120071 ) . is that the forward 
scattering peak results mainly from the diffraction of the incident light. The 
strength of the diffraction peak and its scattering angle dependence appears 
to depend mainly on the size of the particles and fairly shape independent 
for projected-surface-are a equivalent convex particles in random orientation 
dMishchenko et all liooih . 



Because our normalization of the measured phase function, F 11 (0)/F 11 (3O°), 
at = 30° is rather arbitrary (we could have chosen a different value of 
for the normalization), we scale the phase function as calculated for the 
spherical palagonite particles to the measured phase function. For this, we use 
the following equation, 



F n (0°) _ ^(0°) F n (3°) 
Fn(30°) ~ 



FM3 ) Fn(30°)' 



(7) 



with i 71 ii(0 o )/-Fn(30 o ) the artificial data point at = 0° and the superscript 
"s" indicating the phase function as calculated for the spherical particles. 

We now have data points available across the full scattering angle range to 
evaluate the integral in Eq. ([H]) and to obtain the normalization constant C. 
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The numerical evaluation of this integral, however, appears to be difficult be- 
cause of the steep slope between = 0° and 6 = 3°, where no measured data 
points are available. Therefore, we have chosen a different method to obtain 
the normalization constant C in Eq. This method is based on the expan- 
sion of the measured phase function Fn(6)/Fn(30°), including the added, 
artificial, data points at = 0° to = 180°, as a function of the scatter- 



ing angle into so-called generalized spherical functions (jGel'fand et al.l. 11 963 



Hovenier and van der Med . ll983l ; lde Rooij and van der Stapl . ll984t IHovenier et al 
20041 ). This method of obtaining expansion coefficients from scattering matrix 



elements is explained in detail in Sect. HI The expansion of F 11 (0)/F 11 (3O°) 
yields expansion coefficients a[ (with < /). The first of these expa nsion co- 



efficients, a?, is equal to constant C in Eq. ([6]) (IHovenier et all 12004 ). Having 
obtained C in this way, and thus ^"(30°) 
measured ratio F 11 (0)/F 11 (3O°) and Eq. (] 



we readily find F^(Q) from the 



Next, given the auxiliary phase function, we derive the synthetic matrix el- 
ements F^ u (0) for i,j — 1 to 4 with the exception of i — j — 1 from the 
measured ratios Fy(0)/ Fn(0), using Eq. ([3]). To obtain also the complete 
scattering angle range for the other scattering matrix element ratios, we ex- 
trapolate the measured Fy(0)/Fn(0) (z, j = 1 to 4 with the exception of 
i — j — 1) towards = 0° and 180°. At these t wo scattering angles, the 
following equalities should hold: (see Display 2.1 in IHovenier et al.l . 120041 ) 



Fi2(0 o )/F n (0 o ) = F u (0°)/F u (0 ) = 0, (8) 

F 2 2(0°)/Fn(0 o ) = F 33 (0°)/F n (0°), (9) 

F 22 (180 o )/F n (180 o ) = -F 33 (180°) / ^(180°), (10) 

F 12 (180 o )/F n (180 o ) = F 34 (180°) / ^(180°) = 0, (11) 

F u (180°)/F u (180°) = 1 - 2F 22 (180 o )/F n (180 o ). (12) 

Following Eq. ©, ratios F 12 (0 o )/F n (0 o ) and F 34 (0°) / F n (0°) are set equal to 
zero. Following Eq. ([9]), we use splines to extrapolate the ratios i 7 22 (0)/i r n(0) 
and F 33 (0)/F n (0) towards = 0°, and we set both F 22 (0 o )/F n (0 o ) and 
-^33(0°)/i 7 'ii(0 o ) equal to the average of the two extrapolated values. Ratio 
i 7 44(0°)/F 11 (0°) is obtained by extrapolating (with splines) the ratio i 7 44(0)/i r n(0) 
from = 3° towards = 0°. 

In the backward scattering direction, the measured scattering matrix element 
ratios Fj j (0) /F^jQ ) appe ar to be smooth functions of (see Fig. @J. We use 



splines flPress et al.l . Il992f ) to extrapolate F 22 (0)/F n (0), and F 33 (0)/F n (0) 
from = 174° to = 180°. Because F 22 (180 o )/F n (180 o ) should be equal 
to -F 33 (180 o )/Fn(180°) (see Eq. (fTD]) ). we set both F 22 (180°) / F n (180°) and 
— F 33 (180°)/Fn(180°) equal to the average of the two extrapolated values. Fol- 
lowing Eq. dnj, we set F 12 (180 o )/Fn(180°) and F 3 4(180 o )/Fn(180°) equal to 
zero, and calculate F 44 (180°) / Fn(180 o ) using Eq. flfi} with F 22 (180°) / 'F n (180°) 
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In the followin g, we will refer to the e lements of our auxiliary scattering matrix 
F au (6) as (see IHovenier et allliooih : 



F au (6) 



ai(9) 6i(6) 
6i(6) a 2 (e) 
a 3 (e) 62(e) 
-62(6) 04,(0) 



(13) 



4 Expansion coefficients 



4-1 Definitions of the expansion coefficients 



For use in numerical radiative transfer algorithms, it is often advantageous 
to expand the elements of a scattering matrix as fun ctions of the scatter- 



ing angle into so-called generalized spherical functions (jGel'fand et al.l. 1 1963 



Hovenier and van der Med . Il983l ; Ide Rooij and van der Stapl . 1 1984 IHovenier et al 



20041 ) . The advantage of using the coefficients of this expansion, the so-called 
expansion coefficients, instead of the elements of a scattering matrix them- 
selves is that it can significantly speed up multiple scattering calculations, 
which is of par t icular importance when p olarization is taken into account 
dHovenier et all . I2004J : Ide Haan et~aD . Il987h . 



We indicate the generalized spherical functions by P^ n (cose) with the in- 
dices m and n equal to +2, +0, -0, or -2, and with I > max{|m|, |n|}. Note 
that generalized spherical function Pq is simply a Legendre polynomial. The 
expansion of the elements of auxiliary scattering matrix F au (e) (see Eq. ffl3l) ) 
into generalized spherical functions is as follows: 



ai(6) 


= £ai p o,o(oose), 

1=0 


(14) 


a 2 (Q) + a 3 (Q) 


00 

= E(«2 + «3)^ 2 (cose), 

1=2 


(15) 


a 2 (6)-a 3 (e) 


00 

= £(«2-a 3 )4-2(cose), 

1=2 


(16) 


a 4 (6) 


00 

= E«4^0,o( C OS@) ; 
1=0 


(17) 


61(e) 


00 

= E/M, 2 (cose), 


(18) 



1=2 
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M©) = E/«, 2 ( COS0 )- 



(19) 



1=2 



Here, a[, a 2 , az 3 , a4, (3{, and f3 2 are the expansion coefficients. For each value of 
integer Z, the expansion coefficients can be derived from the auxiliary scattering 
matrix elements using the definiti ons of the generalized spherical functions and 



their orthogonality relations (see iHovenier et all 120041 ) . 



A similar expansion in generalized spherical functions can be made for any 
scattering matrix of the form given by Eq. ([2]). The coeffici ent a® is always 



equal to the average of the one-one element over all directions (IHovenier et al. 



2004 ). So for the auxiliary phase function ai(0) we have, according to Eq. ([5]) 
1 and for the measured phase function, F 11 (9)/F 11 (30°) we have a® = 



a; 



C, as mentioned in Sect. 13.41 



4-2 The Singular Value Decomposition method 



To derive the expansion coefficients of the measured phase function, incuding 
the points added at and 180 degrees, and of all six elements of the auxiliary 
scattering matrix we write each of the Eqs. (|T^I) - (|T9l) into the following general 
form: 

m 

y(e) = £7^(©), (20) 

l=n 



where y(Q) represents the value of the measured phase function or an auxiliary 
scattering matrix element at scattering angle O, or, in the case of a 2 and 03, re- 
spectively their sum, as in Eq. ffTBT) . or difference, as in Eq. ffToT) . The functions 
X'(0) in Eq. (I20I) are the basis f unctions, for wh i ch we choose the appropriat e 



generalized spherical functions (IHovenier et all I2004J ; Ide Haan et all 119871 ). 



The parameters 7' in Eq. (1201) represent the expansion coefficients, or, in the 
case of a l 2 and a l 3 , respectively their sum (Eq. (flBT) ) or difference (Eq. (fTB]) ). 
Furthermore in Eq. (1201) . n equals or 2, depending on the auxiliary scattering 
matrix element under consideration (see Eqs. (TT4|) — (TT9]) ) . and, although the- 
oretically m equals 00 (see Eqs. (}T^- (fT9l) ). in practice m is restricted to the 
number of scattering angles at which values of the auxiliary scattering matrix 
elements y(Q) are available. 

For a linear model such as that represented by Eq. (|2"0"|) . the merit function x 2 
is generally defined as: 



E 

i=i 



ym-ET= n YX l (e 



(21) 



17 



where k is the number of available data points and Ui is the error associ- 
ated wit h data point y(Qj) . We use the Singular Value Decomposition (SVD) 
method (IPress et all 119921 ) to solve Eq. (121 j) for the expansion coefficients 7', 
because this method is only slightly susceptible to roundoff errors and provides 
a solution that is the best approximation in the least-squares sense, both for 
overdetermined systems (in which the number of data points is larger than 
the required number of expansion coefficients) and underdetermined systems 
(in which the number of data points is smaller than the required number of 
expansion coefficients). 



To test the robustness and the quality of the fit method based on the SVD 
method, we applied it to scattering matrix elements that we calculated for 
the hypothetical, spherical pa lagonite particles, and that are sh own in Fig. H] 
(note that our Mie-algorithm (jde Rooij and van der Stapl . 119841 ) provides ma- 
trix elements normalized according to Eq. (jSJ), although in Fig. HI the nor- 
malization of the elements has been adapted to correspond to that of the 
measurements). For the test we compare the matrix elements calculated with 
our Mie-algorithm with the matrix elements obtained using the expansion co- 
efficients derived with the SVD method and Eqs. (|14p - (|19p . For the relative 
errors in the matrix elements we adopt the values of the experimental errors. 



We tested two aspects of the application of the SVD method. First, we ap- 
plied the method to matrix elements calculated at the same set of scattering 
angles as the measured ratios of matrix elements, i.e. having a typical angular 
resolution of 5° and scattering angles ranging from 3° to 174°. Comparing the 
matrix elements obtained using the expansion coefficients that were derived 
with the SVD method with the directly calculated matrix elements, we found 
that the relatively coarse angular sampling and the lack of data points below 
= 3° gave rise to strong oscillations in the matrix elements that were cal- 
culated from the derived expansion coefficients. In addition, with the derived 
expansion coefficients, we could not reproduce the strong forward scattering 
peak in the phase function. The lack of data points above = 174° appeared 
to be less of a problem, probably because of the smoothness of the matrix 
elements at those scattering angles. 

Second, we applied the SVD method to matrix elements calculated at an an- 
gular resolution of 1° and covering the full scattering angle range, i.e. from 0° 
to 180°. The matrix elements obtained using the hence derived expansion co- 
efficients coincided within the numerical precision with the directly calculated 
matrix elements. From this we conclude that our implementation of the SVD 
method is reliable, but that we have to apply it to the whole scattering angle 
range, and with a relatively high angular resolution. 

Finally, we found that averaging two sets of derived expansion coefficients, 
one of which has one coefficient more than the other, removes most of the 
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left-over oscillations in the scattering matrix element that is calculated with 
the coefficients. 



4-3 The derived expansion coefficients 



For deriving the expansion coefficients of the scattering matrix elements of 
the Martian analogue palagonite particles, we use the elements of the auxil- 
iary scattering matrix F au (G) (see Sect. 13.41) . because they cover the whole 
scattering angle range and are normalized according to Eq. 

We increase the angular sampling of the auxiliary scattering matrix by adding 
artificial data points by spline interpolation to the dataset between = 3° 
and 174°. Starting with 44 measured data points per matrix element, adding 
the artificial data points results in 220 data points for each of the auxiliary 
scattering matrix elements. This amount of artificial data points includes the 
values at the forward and backward scattering angles at respectively = 0° 
and 6 = 180° as described in Sect. 13.41 The 220 datapoints proved to be 
required to be able to follow the steep slopes of ai(0), a 2 (9), a 3 (Q) and 
04(0) between 6 = 0° and = 3°, where no artifical datapoints were added, 
without introducing unwanted oscillations. As 220 datapoints are needed to 
obtain the auxilliary phase function aj.(G), it proved to be practical to extend 
also i 7 i2(0)/-^ii(0) and F 34: (Q) / Fu(Q) to the same amount of 220 datapoints, 
to be able to straightforwardly apply Eq. [3] to obtain b\ (0) and b 2 (Q). 

The optimal number of expansion coefficients for each of the elements was ob- 
tained in an iterative process: unrealistic oscillations at large scattering angles 
are suppressed when using a smaller number of expansion coefficients, while 
the fit of the steep phase function near 0° is improved when using a larger 
number of expansion coefficients. Applying our SVD method to the auxiliary 
scattering matrix elements, leaves us with respectively 185 expansion coeffi- 
cients for ai(0), 02(0), 03(0), 04(0), 130 expansion coefficients for b± (0) and 
46 expansion coefficients for 62(0)- Figure [5] shows the expansion coefficients 
a\, a l 2 , a l 3 , a l 4 , (3[, and f3 l 2 , derived from the auxiliary scattering matrix of the 
Martian analogue palagonite particles. The expansion coefficients are plotted 
with error bars that originate from the experimental errors in the measure- 
ments. An electronic table of the expansion c oefficients wil l be availabl e from 
the Amsterdam Light Scattering Databas^jO (jVolten et al. , 2005 . 2006 ). 



Website: http: / / www.astro.uva.nl/scatter 
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4-4 The synthetic scattering matrix 



Employing Eqs. ffT^|) - ffTU|) with the expansion coefficients presented in Sect.[ 
we can now calculate, at an arbi trary angular resolu tion, the so-called syn- 
thetic scattering matrix (see also iMunoz et all 12004] ) which covers the com- 
plete scattering range, i.e. from 0° to 180°, and which is normalized according 
to Eq. ([5]). Figure [6] shows the calculated synthetic scattering matrix elements. 
The elements of the synthetic scattering matrix are also listed in Tabled] at an 
angular resolution of 1 to 5 degrees. An electronic table of the synthetic scat- 
tering matri x elements wi l l be availab le from the Amsterdam Light Scattering 
Database 3 (IVolten et all 120051 . 120061 ). As a check we used the synthetic scat- 
tering matrix to compute the same ratios of elements as have been measured. 
We found the differences to lie within the ranges of experimental uncertanties 
or very nearly so. 



5 Summary and discussion 



We present measured ratios of elements of the scatt ering matrix of i rregu- 



larly shaped Mar tian analogue palagonite particles (IRoush and Belli . 11995 



Banin et all 119971 ) as functions of the scattering angle O (3° < < 174°) at 
a wavelength of 632.8 nm. Our measured ratios of scattering matrix elements 
differ strongly from those calculated for homogeneous, spherical particles with 
the same size and refractive index. In particular, the measured phase function 
(ratio -Fii(0)/Fn(3O°)) shows a very strong (almost three orders of magni- 
tude) forward scattering peak and a smooth drop-off towards the largest scat- 
tering angles, where the phase function of the spherical particles shows much 
more angular features especially in the backward scattering direction. Clearly, 
using scattering matrix elements that have been calculated for homogeneous, 
spherical particles when irregularly shaped particles are to be expected in 
radiative transfer calculations for e.g. the interpretation of remote-sensing ob- 
servations can thus lead to errors in retrieved dust properties (microphysical 
parameters and/or optical thicknesses). 

To facilitate the use of our measurements in radiative transfer calculations 
for e.g. Mars, we have first constructed an auxiliary scattering matrix from 
the measured scattering matrix. This auxiliary scattering matrix covers the 
whole scattering angle range (i.e. from 0° to 180°), and its elements have 
been normalized such that the average of the phase function over all scat- 
tering directions equals unity. The value of the phase function at = 0° 
has been computed from the phase function calculated with Mie-theory for 
homogeneous, spherical particles with the same size and composition. The 
normalization of the auxiliary phase function, and hence the normalization 
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of the other elements too, is obtained by applying a Singular Value Decom- 
position (SVD) method to fit an expansion in generalized spherical functions 
( jGeFfand et all Il963l ; iHovenier et all 120041 ) to the measured phase function, 
including artificial data points. The first expansion coefficient yields the re- 
quired normalization constant. After the normalization of the auxiliary phase 
function, the SVD method is applied to the auxiliary scattering matrix and 
its expansion coefficients are obtained. 



With the expansion coefficients, a synthetic scattering matrix is computed for 
the complete scattering range. It is normalized so that the average of its one- 
one element over all directions equals unity. The synthetic scattering matrix 
elements can also straightforwardly be used in radiative transfer calculations. 
The need to include all scattering matrix elements, instead of only the phase 
function, is obvious for the interpretation of polarization observations. How- 
ever, even for flux calculations, all scattering matrix elements should be used, 
because ignoring polarization, i.e. using only the phase functio n, in mufti 
pie sc attering calculations induces errors in calculated fluxes (e.g. lLacis et al. 



19981 ). The use of only the phase function should be limited to single scattering 
calculations for unpolarized incident light. 



Figure [7] shows a comparison between our sy n thetic phase function and two 
phase functions presented by iTomasko et al.l (119991 ) as derived from diffuse 
skylight observations of the Imager for Mars Pathfinder. Ea ch of the phase 



functi ons has its own normalization. The phase functions of ITomasko et al. 



( 119991 ) show the same general angular behavior as our synthetic phase function: 
a strong forward scattering peak and a smooth drop-off towards the largest 
scattering angles. The forward scattering peak of our phase function appear s 
to be stronger than the peaks of the phase functions of ITomasko et al.l ( 19991 ). 
This can easily be due to the size difference of the particles. The dust particles 
in our sample haye an eff ective radius that is a factor of 2 to 3 larger than that 
of Tomasko et al. (Il999h (4.46 /im versus 1.6 fim ± 0.15 /im). The slopes of 
the phase functions in the backward scattering direction are very similar, and 
appear to be typical for (terrestrial) i rregularly shaped mineral particles with 



mode rate refractive indices (see e.g. IVolten et al. 



2001; Munoz et al.. 2000 



20011 ). This smooth slope can not be easily anticipated from Mie-theory. 



The expansion coefficients, the synthetic scattering matrix elements, and the 
measured ratios of elements of the scattering matrix of the Martian analogue 
palagonite particles will all be available f rom the Am sterdam Light Scattering 



Database (for details, see IVolten et al.l . 12005k 120061 ). The Amsterdam Light 



Scattering Database contains a collection of measured scattering matrix el- 
ement ratios, including information on particle sizes and their composition, 
for various types of irregularly shaped particles. The SVD method presented 
in this article can straightforwardly be applied to measured scattering matrix 
elements of particles other than the Martian analogue palagonite particles. 
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Fig. 1. Scanning electron microscope (SEM) image of the Martian analogue palag- 
onite particles used in this study. Note that SEM images generally give a good 
indication of the typical shapes of particles, but not necessarily of their sizes. 
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Fig. 2. Measured size distributions of the Martian analogue palagonite particles. 
The lines refer to the normalized number distribution (plusses), the normalized 
projected-surface-area distribution (triangles), and the normalized volume distri- 
bution (crosses) as functions of logr, with r the r a dius (in microns) of a project- 
ed-surface-area-equivalent sphere (see IVolten et all I2005I . for the definitions of the 
various distributions). Approximating the normalized number distribution with a 
log normal distribution function yields r e g- = 4.46 and v e g = 7.29. 
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Fig. 3. Sketch of the experimental set-up. The light source (a HeNe laser) emits a 
beam of quasi-monochromatic light, which passes through a polarizer and a modu- 
lator (not shown). The light is scattered by the ensemble of palagonite particles. A 
detector measures the flux that is scattered over a scattering angle 0. Measurements 
can be performed for 3° < 6 < 174°. 
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Fig. 4. Ratios of scattering matrix elements as functions of the scattering angle. 
Crosses show the ratios of elements measured, at A = 632.8 nm, for the Martian 
analogue palagonite particles, with the vertical bars indicating the experimental 
errors. Solid lines show the ratios of elements calculated using Mie-theory for ho- 
mogeneous, spherical particles that are distributed in size as shown in Fig. 
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Fig. 5. The expansion coefficients a\, a 2 , a 3 , a 4 , [3\ and f3 2 as functions of integer I, 
with their absolute errors, as derived from the auxiliary scattering matrix elements 
ai(0), a,2(@), 03(6), 04(0), &i(0), and 62(0), respectively. 
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Fig. 6. The synthetic scattering matrix elements a±, a,2, 03, 04, 61, and 62 as functions 
of the scattering angle 0. 
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)K: Mars pathfinder imager results @ 444 nm 
(estimated from figure 21, Tomasko et al. (1999)) 

- : Martian analogue palagonite 




30 



I i i I i i I 

60 90 120 

Scattering angle © (degrees) 



150 



180 



Fig. 7. The synthetic phase function of the M artian analogue palag onite parti- 
cles compared to the phase function derived by iTomasko et al.l ( 19990 from Mars 
Pathfinder results. Each of the curves has its own normalization. 
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Table 1 

The synthetic scattering m atrix elements of the Martian analogue palagonite parti- 
cles (IRoush and Bellll995h . as functions of the scattering angle (in degrees) . The 
scattering angles are identical to those used in the measurements, extended with 
= 0°, 1°, 2° and from 175° to 180° degrees. An electronic ver sion of this table wi ll 



be available from the Amsterdam Light Scattering Database (jVolten et al 
at http:/ /www. astro. uva.nl/scatter. 



20051 ). 
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Table 1 
Continued. 
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